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Abstract 

Describing by a Hubbard type of model a thin armchair graphene ribbon in the armchair hexagon 
>>■ 
^ \ chain limit, one shows in exact terms, that even if the system does not have flat bands at all, at 
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low concentration a mesoscopic sample can have ferromagnetic ground state, being metallic in the 
same time. The mechanism is connected to a common effect of correlations and confinement. 



I. INTRODUCTION 

Carbon-based nanoscale structures holding hexagonal repeating units and different type 
of boundaries attract considerable interest by unprecedented application possibilities in the 
design of organic nanodevicesii^. From these, nanoribbons are of particular interest^, being 
studied especially because of the emergence possibilities of itinerant ferromagnetisnt^"-, this 
subject being also driven by the aim to produce magnetic behavior in organic materials not 
containing magnetic elements^. 

In these systems the shape of the edges is either of zig-zag or armchair type^. The former 
case is usually associated to localized edge states^ which, due to their high degeneracy, 
can lead to flat-band ferromagnetisnt^ (in most cases a spin-polarized many particle, but 
localized ground state). The localized edge states for zig-zag edges on thick samples have 
been observed experimentally by scanning tunneling microscopy—, but it is also known that 
in the thin ribbon case (the chain limit) the localized ferromagnetic nature is no more present 
in exact terms, and the spin polarized state becomes conducting in the low concentration 
limit for mesoscopic samples^i. Nevertheless, the emergence possibility of the edge localized 
states for zig-zag boundaries directed the attention to the armchair edges in the search for 
itinerant ferromagnetism in these structures. 

For armchair carbon-based nanoribbons the bulk states and the localized end states are 
considered entangled by the short-range Coulomb repulsion^. In these conditions, at the 
level of theoretical predictions, itinerant carriers in dispersive bands are supposed to medi- 
ate via exchange coupling the interaction among the local magnetic moments present in flat 
bands^*^. Different approximate techniques are used for the description as: effective field 
theories combined with variational wave function approach^, or analytical weak-coupling 
analysis combined with numerical density matrix renormalization-group method and first- 
principles calculations- . As a general observation one notes that the deduced results attract 
the attention to the importance of the Coulomb interaction and correlation effects in provid- 
ing the physical properties of the systems under study^*^*^. This information is important 
also in a broader context, since for example, potential nanoribbon applications in building 
up organic transistors or spin qubits^i^^ do not take into account correlation effects at all. 

On the theoretical side the studies directed to the explanation of ferromagnetism are in a 
stage present at other organic systems as well, for example pentagon chains^, where the initial 



attempts based on flat band ferromagnetisnt^ encounter the difficulties of this model. One of 
these is the fact that a perfectly flat band is difficult to obtain. The second difficulty is related 
to the observation that the connectivity condition needed for the flat band ferromagnetism, is 
often not satisfied in experimentally manageable situations'^, and this case is encountered for 
armchair nanoribbons as well^. The corrected theory in the light of this new input is usually 
a some kind of developed flat band theory which does not need the in-flat-band connectivity 
condition. On this line, for example in the case of the pentagon chains, a model applied in 
describing rare earth compounds has been suggested to work— using strong hybridization 
effectsi^, while in the case of the armchair nanoribbons constructed from hexagon cells, 
electrons from dispersive bands are conjectured to interconnect the moments created inside 
the flat bands^i^. In fact, this last possibility is also known in the literature. In D = 1 one 
dimensional case it originates from Ref.— , whose extended version to D > 1 can be find in 
Rem. 

All the previously mentioned results relating ferromagnetism in organic nanoribbons sug- 
gest that spin-polarized states in such systems are intimately connected to bare flat bands. 
This suggestion focuses the research directed to the understanding of ferromagnetism (and 
this is the case also for graphene structures) exclusively on flat band emergence possibilities, 
in condition when it is known that systems built up from hexagonal cells containing several 
sublattices, usually do not contain flat band9ii»22^ and flat bands are possible to appear only 
in special conditions, for example as consequences of edge states. Given by the importance 
of the application possibilities of such organic systems in various fields of nanotechnology, 
this state of facts naturally leads to a main question: are indeed flat bands necessary for 
the emergence of ferromagnetism in structures built up from repeating hexagonal units ? 
Guided by this question, with the aim to provide relevant information and advancement in 
this field, one shows below in exact terms, that the answer to this question is negative. 

In order to do that, one chooses the simplest possible armchair nanoribbon, an armchair 
hexagon chain (see Fig.l) which besides the fact that represents a real existing structure 
(i.e. polyphenanthrene) , represents as well an armchair nanoribbon in the extreme thin 
limit. One uses for it a realistic Hubbard type of model which takes into account as well 
the non-zero next nearest neighbor hopping amplitudes^S^ and different on-site potentials 
on different type of sites. Even if the system is not integrable, using a special technique 
one deduces exact ground states for it. These, in the small concentration limit turn out 



to be itinerant and ferromagnetic in conditions in which flat bands are completely missing 
from the system. One shows that the spin polarization emerging for mesoscopic samples is 
a consequence of the common effect of correlations and conflnement. Since increasing the 
ribbon width there is not present a physical reason to cancel this behavior, we expect similar 
effect to occur for thicker ribbons as well. 

The method one uses is related to positive semidefinite operator properties. Based on 
these, one transforms the Hamiltonian of the system [H) in positive semideflnite form, the 
ground state being provided by the most general wave vector which gives the minimum 
possible eigenvalue (i.e. zero) of the positive semideflnite operators building up H. The 
technique itself works well even in unexpected circumstances from the traditional point of 
view of exact solutions as: three dimensions^i, disordered and interacting systems in two 
dimensions'!^ stripes, checkerboards and droplets in two dimensions'^, delocalization effect of 
the Hubbard interaction in two dimensions'^, ferroelectric systems'^, superconductors'^''^, 
or non-integrable chain structures even under the action of external magnetic flelds'^''^. 
Different procedures of the method can be flnd in Refs.— ~— , and extreme details regarding 
its application in the case of the chain structures have been presented in Refs.— "2^. 

One notes that the flrst exact results for the armchair hexagon chain in the interacting 
case are contained in this paper. In order to deduce these, a special technique able to 
handle the emerging extended operators in the presented case has been developed, used, 
and described. 

The remaining part of the paper is organized as follows. Section II. presents in details 
the studied system. Sect. III. describes the transformation of the Hamiltonian in positive 
semideflnite form. Sect. IV. deduces the ground state wave function and presents the phys- 
ical properties of the obtained ferromagnetic phase. Sect. V. containing the summary and 
conclusions closes the presentation, and the Appendices A-C contain mathematical details. 

II. THE SYSTEM 

The studied system is an armchair hexagon chain whose cell constructed at the site j is 
presented in Fig.l. One has 8 in-cell positions at the sites j + r,y, u = 1, 2, ..., 8, a represents 
the Bravais vector of the chain, and for mathematical convenience one has ri = 0. One 
further underlines that the u index denotes as well eight different sublattices S^ present in 
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FIG. 1: The armchair hexagonal chain with cell constructed at the site j (thick line). The in-cell 
positions are j + rj^, where z^ = 1, 2, ..8 is denoting as well eight different sublattices. a is the Bravais 
vector, and one has ri = 0. The numbers inside the cell represent the v index of the in-cell sites. 



the system. 

The Hamiltonian of the system can be written as 



H = To + Ti+T2 + Hu, 



(1) 



where the kinetic part Hq = Tq + Ti + T2 contains the contributions of the one-particle local 
potentials (Tq), the nearest neighbor (Ti) and the next nearest neighbor (T2) hopping terms, 
while Hjj represents the Hubbard interaction. One has 
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In (111) q^ creates an electron with spin a at the site i, h\^a is the particle number operator, 
Nc represents the number of cells, and during the summation over n, j = i + raa is considered, 
i being an arbitrary site of the sublattice Sy=\. The parameters of H describing a realistic 
system are presented in Fig.2. The nearest neighbor hoppings are ti (touching bonds between 
hexagons), t^ (horizontal hopping on the armchair), ts (the oblique hopping on the armchair), 
and ^4 (the hopping on the external armchair bond). The next nearest neighbor hoppings 
are t\ (bonds with only one end on touching points between hexagons), and t'2 (bonds with 
both ends on touching points between hexagons). Furthermore, one notes that the one- 
particle on-site potentials (ea) and the on site Coulomb repulsions \Ja are denoted by the 
index a = on the touching points between hexagons, and by the index a = 1 on other 
(external) sites. For mathematical simplicity one considers f/o = f/i = f/ > 0, and periodic 
boundary conditions are taken into account. 
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FIG. 2: The Hamiltonian parameters (see text following dS])). 

III. THE TRANSFORMATION OF THE HAMILTONIAN IN A POSITIVE 
SEMIDEFINITE FORM 

A. The defined block operators 

In order to transform the Hamiltonian ( TT|2i) in a positive semidefinite form one defines 
eight block operators ^p,j,(T, P = 1, 2, ..., 8 for each site j of the sublattice z/ = 1, as follows 
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(3) 



As seen from ([3]) the used block operators Ap j g. are linear combinations of the starting 
fermionic operators Ci^o- acting on the sites of the block p defined at the site j (see FigslMTOl). 
The numerical prefactors ai,bi,di, ...,ki are unknown at the moment, and will be deduced 
further on. 




FIG. 3: The block p = 1 depicted in the cell defined at the site j providing the block operator 
^1 j^o-- The thick line shows the block, while the numbers are indicating the i' index of the sites 
whose position is j + r^,. 




FIG. 4: The block p = 2 depicted in the cell defined at the site j providing the block operator 
^2,j, cr- The thick line shows the block, while the numbers are indicating the i' index of the sites 
whose position is j + r^,. 

B. The positive semidefinite form of the Hamiltonian 



Based on the block operators defined in ([3]), the starting Hamiltonian presented in (ITP|) 
is transformed (i.e. is rewritten in exact terms) in the following positive semidefinite form 
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FIG. 5: The block p = 3 depicted in the cell defined at the site j providing the block operator 
^3 j^o-- The thick line shows the block, while the numbers are indicating the i' index of the sites 
connected exclusively to a single cell. 
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FIG. 6: The block p = 4 depicted in the cell defined at the site j providing the block operator 
^4 jo-. The thick line shows the block, while the numbers are indicating the i' index of the sites 
connected exclusively to a single cell. 



where A^ represents the operator of the total number of electrons and K is a constant. The 
numerical prefactors in ([3]) and the constant K are expressed by i) effectuating the products 
Ap ; ^Ap j o- and the prescribed sums in the first term of (jl]), and ii) equating the obtained 
result term by term to the expression of the starting Hamiltonian (TTJI^]) . The obtained 
equations connecting the parameters of the starting Hamiltonian to the parameters of the 
block operators and the constant K in the transcribed Hamiltonian (jl]) are the following 
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FIG. 7: The block p = 5 depicted in the cell defined at the site j providing the block operator 
^5 jo-- The thick line shows the block, while the numbers are indicating the i' index of the sites 
whose position is j + r^. 
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FIG. 8: The block p = 6 depicted in the cell defined at the site j providing the block operator 
^ej^o-- The thick line shows the block, the number four indicates the z/ = 4 site. Other block sites 
(placed in the following cell) are indicated by their position j + rjy + a. 
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FIG. 9: The block p = 7 depicted in the cell defined at the site j providing the block operator 
^7j^o-- The thick line shows the block, the number eight indicates the u = 8 site. Other block sites 
(placed in the following cell) are indicated by their position j + r,^ + a. 
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FIG. 10: The block p = 8 depicted in the cell defined at the site j providing the block operator 
yl8,j,o-. The thick line shows the block, while the numbers are indicating the v index of the sites 
whose position is j + Fj^. 
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ones: The hopping terms become expressed as 



^j+r5,j+r4 = ^1 = ^5^4 + dld^ + f^ fi + k^k^, tj+nj+ro = ti = ^i^e + go96 + ^0^6 + d^de, 

^j+rej+rs = ^2 = ^6^5 + ^6^5 + ^6^5? '^j+ri+a,j+r4 = h = d^d^ + /q /4 + 5'o5'4, 

^j+r2,j+ri = h = OaOi + 5'25'0, ^j+r4,j+r3 = ^3 = «4'^3 + /4/3, 

ij+re+aj+rg = ^3 = CgCs + /ig^S, ^j+ryj+rs = ^3 = 6765 + kjk^, 
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'■j+rsJ+ra ~ n ~ Js/S) ^j+r4J+r2 = n ~ '24'3.25 '^j+r2J+r6 = n ~ 5'25'6) 

''j+rgj+ri+a = ^-l = '^s'^Oj '-j+rsj+rg ~ '^l ~ ^5^85 ''j+r4,j+r7 ~ ''I ~ "^4"^7) 

'"j+rrj+re ~ n ~ kjkQ, tj+rrj+rB+a = n ~ ^7^65 i^j+n+aj+rs = n ~ /o/S) 

^j+r2+a,j+r4 = i^l = 5'25'4, ^j+rs+aj+rg = ^1 = ^5^8, ^j+r6+a,j+r4 = ^2 = d^di + 5'g5'4, 

^j+rij+rs = ^2 = ^1^5 + "-Q^'S) ^j+r6,j+r4 = ^2 = "6"4 + H^45 



t 



'j+ri+aj+rs 



^2 — '^O'^S + /o/s) 



ij+re+aj+rg — — (ig(i5 + 6565, tj+r4J+ri — — 0401 + 6461, 

while the on-site potentials are given by 
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flip + i&ip + Mor + i/or + i^or + i^or + k, 

|2 I Ii,l2 I |j|2 I lxl2 I I |2 I 1; |2 I r^ 

04] + I&4I + |d4|^ + 1/4!^ + ig^f + {kil + K, 

&5I' + 141' + lesP + I/5I' + l/isP + \k,\' + K, 

&6p + l^r + leer + l^er + l^eT + \h\^ + K, 
er\^ + \k7\^ + K, 

e8p + \hs\'^ + K. 



(5) 



(6) 



For clarity, in ([5]) one specifies also the hopping matrix elements ti^.ia, as the component of 
^ii,i2Ci 0-^12,0- kinetic energy term which leads to the respective equation. Again for clarity, 
at the level of the on-site potentials in (jHD, in the index of ei, the site i where the analyzed 
potential is considered is precisely indicated. 

One notes that Eqs. (!5]|6|) are the matching conditions of the solution since the transfor- 
mation from (IT1[2|) to (jlj) is valid only if Eqs. (I5]|6l) are satisfied. 
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C. The explicit expression of the operators ^pj,cr 

The exphcit expression of the transformed Hamiltonian in @, hence also the exphcit 
form of the Apj o- operators, can be find by solving the matching conditions f l5|6|) . These 
are representing 36 coupled nonlinear complex algebraic equations containing 33 unknown 
variables, namely the numerical prefactors Oj, 6j, (ij, ..., fcj in ([3]), (in total 32 variables), and 
the constant K in (^. Based on the solution deduced in the Appendix A, see (1A13I) . the 
introduced block operators have the expression 

t' 

^l,j,cr = Vt4 [ (Cj+r2,o- + Cj+rg.cr) + —{C},a + Cj+r4,o-) ]) 

A I — r / 2 1 ~^ -^^2 \ / -^ -^ \ 1 / -^ ^ \ 1 

^2j,<T = V^4 [ (t7 17r )(Cj+r5,cr — Cj+rg.CTJ + — (Cj,^ — Ci^r^.a) J, 

I — r / 2 1 I 1 \ / -^ ■^ \ 1/-^ ■^ \ 1 

^3j,cr = V^4 [ (77 TT: )(Cj+a,cr — Cj+r4,o-) + — (Cj+r5,cr " Cj+re+a.cr) J, 

t' 

^4j,o- = yti [ (Cj+r7,o- + Cj+rs.cr) + " (Cj+rs.o" + Cj+re+a.cr) J, (7) 
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A / / 1 "•" 2 r / -^ -^ N 3 '~ V 3 — 1/ 



-j+a,(T J 5 



4 _ / / ^1+2^2... . 2t3 „ (^3-^'l). 



A / / 1 "•" 2 r / " " N 3 " V 3 — 1/ 

1 vS J- ' 2 1 
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The parameter space domain where this solution of the matching conditions is valid is not 
severely restricted, being given by ^4 > 0, t^t^/^ti + 21^) > 0, and flAlip . see Appendix A. 



The constant K present in the transformed Hamiltonian from @ is given in this case by 

K = e,-U--^^^. (8) 

{ti + 2Q ^ ' 

One notes that the Apj^o- operators present in ([7]) are Fermi operators, since A ■ o-^I ; a ~ 
^p,j,(T^p,j,(T = is satisfied. But these operators are not canonical Fermi operators, since 
{^J,-'^J',j',a'} 7^ '^py'^JJ'^-,-' holds. 
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IV. THE DEDUCTION OF THE GROUND STATE WAVE FUNCTION 

A. The deduction method 

Based on fITfS]) now one knows the exphcit expression of the transformed Hamiltonian 
(jl]), hence one starts the deduction of the ground state wave function for it, which is written 
in the form 

\^,) = [ll^im. (9) 

ai=I 

In ([H]), |0) represents the bare vacuum, Bj^ are hnear combinations of the starting canonical 
Fermi creation operators q^, a represents an index which labels different i?^ operators, and 
/ represents the domain from which the a indices must be choosen in order to build up the 
ground state wave function. 

The strategy in deducing the Bj^ operators is based on the positive semidefinite structure 
of H from (^, namely on the fact that H' = H — KN = Ha + Hjj is a positive semidefinite 
operator, where Ha = Y.a Y.l=i EjI'i ^A,p,j,a and i^A,p,j,a = ^pj,a^p,j,a- Consequently, |^g) 
is the most general wave vector satisfying the requirement 

{HA + Hu)\^g) = 0, (10) 

which, because both Ha and Hu are independently positive semidefinite, is satisfied only if 
both HAl'^g) = 0, and Hul'^g) = hold. But since the kernel ker{0), of an operator O, is 
a Hilbert subspace spanned by all vectors |0) satisfying O|0) = 0, it results that in order to 
deduce |\E'g) we must deduce the most general wave vectors placed inside 

ker{H') = ker{HA) n ker{Hu). (11) 

Concentrating first on ker{HA), one observes that since HA,p,},a are all positive semidefi- 
nite operators for all values of all indices, it results that 

ker{HA) = ker{HA,p=i,i,,^) n A;er(i^A,p=i,ji,i) n ker{HA,p=2,h,^) ^ ker{HA,p=2,ii,l) ^ •••• 

... n /cer(i^A,p=8,j^„t) n ker{HA,p=8j^^,i), (12) 

where, in the right side of the equality, all indices p, j, o" occur. 

Furthermore, one knows that for fixed p = p',j = j', cr = cr' indices, ker{HA,p',y,a') is 
spanned by vectors of the form l"^ a,p' j' y) = Ylaei , ., , ^a\^) where lA,p',y,a' collects all a 
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indices for which the anti-commutation relation {Ap/ j' g-/, S^} = holds^^. Indeed in this 

case HA,p'j',a'\'^A,p',y,a') = '^,',y,a'^p',yyllaei , ., ,^11^) = is satisfied, since based on 
the prescribed anti-commutation relation Ap/ j/ o-/_B^|0) = — i?^Ap/ j/ o-/|0) = 0, which given 
by y4p/j/_o-/|0) = 0, is indeed a true relation. 

It results that based on (TT^ . the vectors contained in ker{HA) can be deduced by deter- 
mining 5^ from the relation 

{ip,j,.,5t} = 0, (13) 

where (IT^ must be satisfied for all values of all indices. One notes that different S^ solutions 
of (fT3|l must be linearly independent. 

After obtaining the B^. operators from flT3l) . the domain / present in ([9]) is such fixed to 
push |\E'g) also in ker{Hu). This last kernel is simply the Hilbert subspace containing all 
wave vectors with zero double occupancy. 

B. The calculation of B'a operators 

1. The I3a expression 

In deducing the Bj^ operators one starts from flT^ and the already deduced expressions 
of the Apjo- block operators presented in ([7]). One considers 

Bi= Yl ^-.i4.>' (14) 

where two kind of unknown variables must be determined, namely i) the domain 'Db^ on 
which the operator S^ is defined, and ii) the numerical prefactors x^ i which specify the 
linear combination on which B^, is constructed. After fixing the B^. expression in flT^ . the 
equation (TT3l) must be written for all Apj o- block operators, obtaining a coupled system of 
16A'"c equations which must be solved. 

In performing this job one observes the following properties, a) The domain T>b^ cannot 
be restricted (e.g. extended over Ni, < Nc cells). The reason for this is that for restricted 
T>B^, we always find at least one Apj o- operator (holding some fixed indices p = p', j = j', o" = 
a') which touches the i?^ operator only in a single point j*. In this case, since only one kind 
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of starting fermionic operator is acting on each site and {cj*,o-,ct } = 1 holds, it results 
that {Apij'y, Bj^} = cannot be satisfied. Consequently, the -B^ operators will be extended 
operators (i.e. the domain P^^ will be extended along the whole system), b) For similar 
reasons it is not possible to have in the expression of Bj^ operators islands of cI, ^, operators 
in the background of c] ^ operators where a' 7^ a. Consequently, for a given B'l in ( 1T^ all 
(Ji indices become fixed, hence ctj = a holds. 

Given by these properties, the a index in (HM specifies in fact a fixed spin index a, and a 
label index i which numbers the possible B^ operators, hence a = {£, a), and the expression 
of the i?^ operators from (IT^ becomes 

4V = E^^4- (15) 

i 

where the i sum extends over all sites of the system (in all sublattices) . 

2. The calculation procedure 

For the effective calculation of the xe^i numerical prefactors of the extended operators as 
bI^ present in flTSj) . special techniques are needed. Such a technique has been developed 
and presented in Ref.— , based on a recursive connection between coefficients in neighboring 
cells described by a non-symmetric matrix R, the prefactors being deduced via the unity 
eigenvalues of i?™, m being an arbitrary positive integer. In order to enlarge and develop 
the techniques applicable in treating extended operators, one presents and uses here another 
alternative method. This is based on deducing the i?j ^ operators for a relatively small chain 
- named test chain - containing m cells, the test chain being closed in a ring by periodic 
boundary conditions. Once these operators are deduced, the obtained result is used in the 
description of a long chain containing M = p x m cells, the long chain being as well treated 
by periodic boundary conditions, p being an arbitrary positive integer. 

In this method first one deduces the linearly independent BJ ^ operators for the test chain. 
This last contains m cells, which are in order numbered by the index Uc = 1,2, ...,m. The 
coefficients xi^i become in this way a;^,nc,i^ where u denotes the lattice site j + rj,, i/ = 1, 2, ..., 8 
inside the Ucth cell defined at the site j of the Si sublattice (see Fig.l). Consequently, in 



16 



this notation, by changing the sum over i in a sum over ric and z/, ( !T5|) becomes 

m 8 



4^,-= $Z SZ^^'^'^j'^+an.+r^,., (16) 



nc.=l !^=1 



where j is an arbitrary site of the sublattice 5*1. Based on fll6p and ([3]), the anti-commutation 
relation flT^ provides the following system of equations for the prefactors of the B\^ oper- 
ators 



aia;£,i,i + a2a;£,i,2 + ctsa^^.i.s + a4a;£,i,4 = 0, 
hxLi,i + &4a;^,i,4 + hxtx^ + &6a;£,i,6 = 0, 
doXe,2,i + d4Xe,i,4: + d^xn^i^^ + dQXi>^2fi = 0, 
e5a;£,i,5 + 66^^,2,6 + 67x^,1,7 + e8a;£,i,8 = 0, 

foX£,2,l + f3X£,l,3 + /4a;£,l,4 + f5X£,l,5 = 0, 
goXe,2A + 92X1,2,2 + fl'4a;£,l,4 + 96X1,2,6 = 0, 

hoXe,2,i + hxe,2,5 + hQXt,2,6 + h8Xe,i,s = 0, 
^43:^,1,4 + hxe,i,5 + kQXe,i,6 + ^7X1,1,7 = 0, 
aiXi,2,i + a2Xi,2,2 + a'3Xe,2,3 + 0-4X1,2,4 = 0, 

blXi,2,l + hiXl,2,4: + hXi>,2,^ + hXl>,2,6 = 0, 
doXe,3,l + '^43^^,2,4 + d5Xi,2,5 + dQXe,3,Q = 0, 

e5Xe,2,5 + eeXi,3,Q + 67X^,2,7 + esXi,2,8 = 0, 

foXe,3,l + f3XL2,3 + f 4X1,2,4 + i^Xi,2,^ = 0, 

9oXe,3,i + 92Xe,3,2 + 94x^,2,4 + 9^x^,3,6 = 0, 
hoXe,3,i + hxe,3,5 + hGXe,3,Q + hsXe,2,s = 0, 
k4Xe,2,4 + hxe,2,5 + kQXi,2,6 + k7Xe,2,7 = 0, 
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O-lXi^m,! + a'2Xe,m,2 + 0-3X1,171,3 + 0,AXe,m,4 = 0, 
blXi^rn,! + b4X£^rn,4 + &53^£,m,5 + &63^£,m,6 = 0, 
C^O^^^,!,! + dAX£^rn,A + C^S^^^.m.S + ^^63^^,1,6 = 0, 

e5Xe,m,5 + e6a;£,i,6 + 67x^,^,7 + esXi^rn,s = 0, 

foXe,l,l + f3Xe,m,3 + f4:Xe,m,4 + f5Xe,m,5 = 0, 

5'oa;£,i,i + 5'2a;£,i,2 + 5'4a;£,m,4 + gaxe^i^e = 0, 
/ioa;^,!,! + hxe,i,5 + hxe^ifi + hsXi^rn,s = 0, 

^4a;£,m,4 + k^Xe^^rn,^ + A;6a;£,m,6 + k-rXg^rn,! = 0, (17) 

where the first eight equations from (fT71) represent (fT3|) written with the operators from ([3]) 
defined in the first cell, the second eight equations from fITTj) represent flT3l) written with 
the operators from ([3]) defined in the second cell, and so on, the last eight equalities from 
flT7|l represent flT^ written with the operators from Q defined in the last mth cell. The 
coefficients ai,...,ki are presented in (lA12p or (jTj). One finds in this manner a system of 
8m linear homogeneous equations for 8m unknown variables, the Xi^ncu prefactors, i being 
fixed. Solutions are present only if the determinant of the system flTTl) is zero. This condition 
specifies different interconnection possibilities between Hq parameters (e.g. different regions 
of the parameter space specified by the index 7) in a generic form 

F,{{t^},{e^}) = 0, (18) 

where solutions are possible to occur. In (ITS]) the sets {t^} and {e^} are representing the 
set of hopping matrix elements, and on-site one-particle potentials, respectively. The study 
of the physical meaning of the condition flTSl) shows that since flat bands are not present in 
the bare band structure (see Appendix B, Eq. (lB7p together with explications related to this 
equation), (fTSi) is not connected to the flat band notion, consequently, the ferromagnetism 
which will be here deduced is not of flat band type. Instead, (IT8II specifles regions in the 
parameter space where the density of states around the minimum system energy value is 
high. 

One given deduced 5^ operator from (TT7|) receives a flxed i index, for example £ = 1. It 
has hence the form 

m 8 
1,0- — / , / ,Xi^n,uC\j^nsi+r„,a i \^^) 



n=l v=l 
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where j is an arbitrary site of the sublattice 5*1. The obtained xi^n,u coefficients in (fT9|) 
usually have the following properties: i) are different in different cells, e.g. for n ^ n' one 
has Xi^nM 7^ Xi^n'M, and ii) rotating the Bl^ operator by n angle around the axis (the line) 
of the chain, the prefactors attached to each site change their value, e.g. 

2^1,n,5 7^ 2;i^„,i, X\^n,& 7^ a^l,n-l,4, a^l,n,7 7^ Xl,n,2, Xi^n,8 7^ Xi^n,3- (20) 

Given by the properties i), ii), new linearly independent BJ^ operators are obtained as 
follows. 

i) One translates all coefficients in flT^ with one cell, the procedure being possible to be 
repeated m — 1 times. One obtains in this manner 

m 8 

n=l v=l 
m 8 



^3,<7 - X. X. ^l'"-2.'^^j4 






m 8 

m,o- ~ / . / ,^'^,'n—m+l,uC^+na+r^,a- (21) 

n=l u=l 

ii) After this step one rotates by vr degree around the axis of the chain all operators 
BI^, i = 1,2, ...,m, obtaining m new linearly independent BJ^ operators for which i = 
m + l,m + 2, ..., 2m holds. The expression of a BI, ^ operator obtained by rotating BJ^ from 
(TTB]) becomes 

m 8 

71c- =1 U=l 
J 

where j is an arbitrary site of Si, and one has (see also fl20l) ) 

/ / / / 

•^l',n,l — -^£,71—1,5: -^l.' ,n,2 — -^f,?!— 1,7) -^l.' .n,Z — ^^,n— 1,8; •^£'.?i,4 — •^^,n,6) 

^£',n,5 ~ 2;£_„^i, %^„_6 ~ Xi,n~l,4, %,n,7 ~ ^i,n,2, %,n,8 = ^£,11,3- (23) 

The linear independence of the new operators obtained via (I21f22j) must be always checked, 
and the presented procedure leads [for one given solution of ( TT7I) ] to maximum 2m linearly 
independent IBJ^ operators for a fixed spin index a. The described procedure than must 
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be repeated for all independent solutions of (fT71) connected to the same parameter space 
domain f lTSj) . Explicit examples for the m = 6 case (twelve hexagons connected in six cells) 
are presented in Appendix C. 

The deduced i?]^ operators on the test chain containing m cells, can be used for the 
description of a chain containing M = p x m cells as well, as stated before, p being an 
arbitrary integer number. For this, a given I3j^ operator (maintaining i and a fixed) is 
periodically repeated p times to cover completely the chain containing M cells (treated with 
periodic boundary conditions). In order to increase the number of linearly independent B^ 
operators relating the same chain containing M cells, one must collect all solutions describing 
the same parameter space region (llSp with different m' values such to reobtain the same 
fixed M via M = p' x m', where both p' and m' are integers. 

C. The ground state wave function 

1. The expression of the ground state 

Collecting together all linearly independent S]^ operators deduced in the previous sub- 
section and denoting their number by A''^, one has for {"^(Nb)) = Ho- Il^Ji -^lo-lO) ^^^ 
property 

^a|*(A^b)) = 0. (24) 

The equality (HM is a direct consequence of ( TT3l) [see also the explications presented just 
above flTSl) ]. Since the kernel ker{HA) is intimately connected to flT^ . it results that in the 
studied conditions, the wave vector \^{Nb)) is unique in satisfying (I2l|) . 

The ground state is obtained from |\1'(A'^b)) based on (TTT]) by taking from it the contri- 
butions which are present also in ker{Hu), i.e. does not contain double occupancy. Since 
all Bj^ operators are extended along a confined space region defined by the studied chain, 
different BJ^ operators holding different i indices intersect each other on several sites of the 
system. Consequently, the double occupancy can be avoided only if one fixes to the same 
value all spin indices of all I3J^ operators. This property remains also true if one decreases 
the number of the i?^ operators present in the wave vector. Hence, the ground state wave 
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function at A^ < Nb becomes 

N 

\^,iN<Ns)) = l[Bl\0), (25) 

e=i 

where a is fixed. The corresponding ground state energy according to (^ is Eg = KN, 

where the K constant is given in (|H]). The ground state nature arises since because of the 

missing double occupancy one has Hul'^giN < Nb)) = 0, while because of f l24|) . i^A|^g(^ < 

Nb)) = 0. Hence indeed {Ha + Hu)\'^g{N < Nb)) = i^'|^g(iV < Nb)) = is satisfied. 

One further note that the described ferromagnetic ground state emerges in the relatively 

low concentration limit. For example, using the explicit solution presented in Appendix C 

in the case of a system with Nc = 36 cells (containing 72 hexagons), the number of electrons 

per cell N/N^ < 0.69 holds, where max{N) = Nb is satisfied. 

2. Physical properties of the ground state 

The deduced ground state fl25l) is a fully polarized ferromagnet. Besides, it is a conducting 
state since the operators building up the ground state wave vector are extended, and 6fi = 
Ij+-I2_ = {Eg{N + l)-Eg{N))-{Eg{N)-Eg{N-l)) = hoMs for A^ < A^^. All itinerant 
electrons being spin-polarized, such systems are potentially applicable in spintronics devices. 

The obtained ferromagnetism is not of flat band type, since as shown in Appendix B, flat 
bands are not present in the system. The spin polarization occurs since extended operators 
conflned in the restricted space region of the chain touch each other, hence must correlate 
their spin index in order to reduce the double occupancy. A such type of ferromagnetism 
is created by the common effect of correlations and conflnement, and was observed in other 
hexagon chains as well^. For the appearance of this phase, besides the conflned geometry, 
the following aspects should be present: a) The operators 5^ building up the ground state 
wave vector must be extended. Such a property is not a priori a chain characteristics, and 
several chain examples are known where it does not occur— i^i^. Regarding this property 
one mentions that always when a such type of behavior was observed for chain structures not 
possessing flat bands, external hoppings (i.e. hopping matrix elements connecting different 
sites from different neighboring cells - see t'l in Fig. 2) were present. Consequently, when 
the hopping terms allow the movement along the whole chain (e.g. cage regions where the 
hopping amplitudes are missing along a chain section perpendicular to the line of the chain 

21 



are not present), this seems to be a sufficient condition preserving the extended nature of 13^ 
operators, b) The density of states at and around the minimum one-particle energy should 
be high. This preserves the presence of the fully saturated ferromagnetic state at arbitrary 
repulsive U in the low concentration limit and finite chain length in the presented case. The 
mechanism however surely works (probably with lower U bound) in the thermodynamic 
limit as well, but the effective proof of this statement remains a challenging problem for 
further investigations. 

One further notes that the described ferromagnetic phase emerges for a sequence of 
different possible interdependences between hopping amplitudes [see the different A = 
(ti — 4:t2)/(ti + 2t'2) values presented after Eq.f lC5|) ]. Since it is known that carbon-based 
nanoribbons are sensitive to uniaxial or shear strains^i, and/or external pressure^, the 
needed interdependences between Hq parameters can be achieved for example by strain 
effects (see also Ref.— ). 

V. SUMMARY AND CONCLUSIONS 

An armchair hexagon chain (i.e. polyphenanthrene) as a representative of a thin arm- 
chair ribbon has been described in exact terms by the use of a technique based on positive 
semidefinite operator properties. The deduced ground state for mesoscopic samples, in the 
small concentration limit was shown to be ferromagnetic and metallic. This appears in con- 
ditions in which flat bands are not present in the bare band structure of the system, the 
spin polarization being provided by a common effect of correlations and confinement. Since 
increasing the ribbon width there are not present physical reasons to cancel in exact terms 
this behavior, at least for a restricted domain of hopping amplitudes we expect similar effect 
to occur for thicker ribbons as well. 
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Appendix A: The deduction of the numerical prefactors of the Apjcr operators 

One starts by using from (|5]) those 14 equations from which every one contains only two 
unknown variables. From these, namely 

ti = a^a2, 14 = 6367, t^ = flgOi, t]^ = J5/3, 

tl — 0402, tl — 9296, tl — '^s'^O, tl — 6568, 

tl = klkj, tl = 6766, tl = /0/35 tl = 9294, 

t'l = k;ke, t'l = hlhs, (Al) 

one finds 



t'l 


a2 = 


t4 

ar 


t'* 
«4 = ^^ as, 

^4 




/5 - .* 

is 


65=, 
^8 


ee = 


t[ 
^ 7768, 

t4 


t4 

67 = -^, 
^8 


^4 = -7, 
i/2 


^6 = -7, 
i/2 


*-.•• 


h = 




A- - 1 


*-| 





(A2) 

Now one uses from (|5]) the 10 equations containing respectively four unknown variables, 
namely 

ts = a*ai + 9^90, t^ = 0^3 + /4V3, ^3 = e^Cg + fc^fcg, 

ts = egCs + hlhs, 4 = &t&5 + ^^^^5, 4 = 65^4 + ^6^4, 

^2 = c^o^5 + /0/5, t'2 = dld4 + 9*94, = 0*01 + 6461, 

= t/*rf5 + e^es, (A3) 

in which, (IA2p is introduced, and provides 



t2% 



{t[y , _ 1 r,/ (O^l U _U ti rt'i 

1 If' |2 + 1+' |2 y-/2 

J _ Ul * ril 1 ^ _ A 4 r|''l| /n J _ ''1 

"0 - ;7r F2 - TjT^J, 04 - 0'5-njT^ [|— [^ - y, de - -TTf^ 
"5 1/3 1 Iti J \92\ t^ar^ 

1 y-' *l/7 |2 1 +'+* 1 y-' U |2 

. _ ^ r,* h |a3| ^ i r, 4^4 , _ ^ r,* HN^t 

/4 - 77 F3 t;; — J, fi'o - -7 F3 - T—i^J, /^6 - T7 Fa -^J, 

J3 '^A 92 P3| "-8 M 

"-7 p8| 
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After this step one remains with 9 unknown variables (03, bi, d^, eg, /s, (72, hs, kj, and K), and 
12 equations of the following form 

h = hlhi + dldi + /s /4 + K^A, ti = blbe + g^ge + h^hQ + d^de, 
t2 = hlh + K^b + ^6^5, h = dldi + /o /4 + 5'o5'4, 

I 12 111,12 11^12 11X12 11 l2ili,l2ir^ 

eo = |a4p + |&4p + M4p + I/4P + |^4p + |fc4p + ^, 

eo = \h? + 141' + lesl' + I/5I' + \h5? + \h? + ^, 

eo = l&eT + |4r + leer + l^el' + M"" + l^el' + K, 

ei = |a2|' + |^2p + i^, ei = lasl' + I/3P + i^, 

ei = |e7|' + |/c7|' + ^, ei = |e8|' + |/i8|' + ^. (A5) 



Using ( JA2JA4I) in ( JA5I) . introducing the notations 



I |2 I J |2 I -P |2 I |2 J. 

lasl =x, Idsl =2;, I/3I =W7, \g2\ = t, 

|e8|' = t', \bi\^ = y, \h8\^ = u, Ikr]'^ = s, (A6) 
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and considering real hopping matrix elements, the remaining system of 12 equations becomes 
of the form 

ei = ^ + t + K, ei = ^ + s + K, 

X V 

ei = X + w + K, ei = V + u + K, 

+'2 ^/2 ^^2 ,,4 , ,/ ^2 ^^2 

1 I 1 I ''1 -^ I ''1 I ,'+ ''1"''n2 I ''4^ /''I ./\2 I TV- 

s t t| tiz u ti t[ s 

,i2 , / 2 , / 2 -| , / 2 -, , / , 

-V I 1 I 1 I 1 I [+' 1 ^2 I J^ /, '-i'-4n2 I TV- 

V w u y us V 

-/ 2 , / 2 , / 2 1 , / 2 -, , / , 

eo = y-\ \ \ ^-1*2 J + 7 1*3 ) + J^ , 

X U W Z W t X 

'- 1 ('- W - 'i' - ^ ('- if ) 4 ('- ^)' 



t4y U t[ t W 14 S V 



4^ 



tiZ W t[ S U 14 t X 

The study of flA7p shows that it provides the same equation for pair of variables. For example 
if rj denotes x or v, one finds the equation 

t t'"^ t'^ 
^2 = 372(^2 -p )( 72 -h) 

+ |^(i3 - ^) + ^(t3 - ^). (A8) 

Similar equations are possible to deduce for the {y,z), {u,w), and (t,s) pairs as well, hence 
one restricts ourselves below to the class of solutions 

X = V, y = z, u = w, t = s. (A9) 
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By using (1A9P in (lA7p . (1A7P reduces to six equations of the form 



ei = '-^ + t + K 



) 

X 



ei = X + w + K, 

./ 2 r, ,; 2 1 ./ 2 J./ 2 , , 
-V I 1 I 1 I /'+' 1 ^2 I ''1 /''3 ''4n2 I TV- 
Co = 2 H \ h - (^2 ) + -- (tt ) + ^' 

X W7 Z W t t\ X 

h = -^ (4 - ^) - ^ (t'2 - ^) + ^ (^ - -) + ^ (^ - -). (AlO) 

One of the possible solutions of flAlOp emerges at x = f = t4 > 0, 2; = y = t'^ /t4 > 0, 
t = u = s = w = 2t3t[/{ti + 2t'2) >0, K = ei-t4- 2Ut\/{ti + 2t^), and requires 

''2\2 






• 2 



+ t;^i^[(|-ir-^ + 2]. (All) 



It has the form 



/' f' 

1 / / 1 

ai = — ;=, (^2 = yti^ c^s = V^4; o-a = ~~/^^ 

V^A V*4 
df) = \/U{Tr TT. )) d4 = -y/U{— — ), rfs = — =, (Iq — 



'^1 2^3 t^ 2^3 yT4 -y/t4 

t' t' 

65 = ^, eg = ^, 67 = v^, eg = v^, (A12) 

V*4 

_ / ti + 2t- _ 



^0 = (^3-^1) 
ho = 



hi + 2t'2 




</w' *^^"'-'' 




26 



The parameter space region V where this solution emerges is fixed by Eq. (lAlip and the 
conditions t^ > 0, t^t'^/{ti + 2t2) > 0. The corresponding Apj,cr operators are provided by 
the following relations 

f' f' 

V''4 V''4 



^1 -^ ^1 ^ / — / ^2 ^1 ~l~ ^^^2 N " I — I 2 ^1 ~l~ ^'^l 



^2,j,a — ~^Cj^ff -=Cj+r4,a + Vt4 (t; ;TT JCj+rg.a — V^4 (t; ;^— " 

VM V''4 f-l '^''3 i^l ^t^S 



2- 



O 



-J+re, cr; 



^2 ^1 "•" ^2 N " / — / 2 ^1 ~r ■^'-'2 \ " 1 -^ 1 



'^j+rSjO- /r—Cj+re+a, cr; 



ti ^t3 r^ /Is 1/14 1/14 

t' t' 

^4j,o- = ^=Cj+r5,(T H ^Cj+r6+a,cr + Vi^4Cj+r7,(T + V t4Cj+r8,o-, (A13) 

A/t4 A/r4 

^ _ , /ti + 2t'2 . / 2t;t3 . /ti + 2t^ . , /ti + 2t'2 . 






/ ti + 2t-2 / ^i + 2t-2 . ,. l h + 2t', 2t;t3 

. /ti + 2t'2 . , , ,, [h + K^ , ,/ /^i + 2t'2 . , / 2t;t3 



Appendix B: The bare band structure 

In order to obtain the bare band structure of the system, one Fourier transforms Hq = 
To + Ti + T2 to k space by using Cj+r,,,^ = (1/v^) Zlk <3xp[-'ik(j + ri,)]c,,_k,f7 (where the sum 
over k represents a sum over Nc cells), obtaining 
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a k 






jkre _|_ st A , +Jk(r5-r4)^ 



^1 = 2^ Z^[^i(cI,k,.C6,k,.e *^^« + c^^^C4,k,.e 

cr k 

-Unfr''' r^, p+*k(r6-r5) I st s^ , +ik(a-r4)N 

+^3(4,k,.ci,k,.e+^'^^^ + cik,.a3,k,.e+^'^(^*-^3) + 4,k,.C5,k,.e+^'^(^-^^) + 4,k,.C8,k,.e+^'^(^«+^-«)) 

+t4(4,k,.C2,k,.e+^'^(^3-r,) ^ .t^^^^.^^^^^^+.k(rs-r,)) ^ ^^^^]^ 

7; _ ^P^Pr+'/At A H-Jkrs , -it A +Jk(r4-r2) , At A +Jk(r5-r3) 

cr k 

_|_At A +Jk(r2-r6) I -t A +Jk(r4-r7) I -t A +ik(r5-r8) 

_|_At A , +ik(r7-r6-a) , -t A +ik(r8-a) , -t A +Jk(a-r3) 

_|_At A +Jk(r2+a-r4) I -t A +ik(r7-r6) , -t A +ik(r5+a-r8)\ 

+4(4,k,.C5,k,.e-'^r'^' + 4,k,.C4,k,.e+^'^(^«-^) + 4,k,.C5,k,.e+^'^(-r'^^) 
+4,k,.C4,k,.e+*^('-«+^-^*)) + H.c.]. (Bl) 

Introducing the notation /ca = ka G (— vr, +7r], b = a/3, where a = |a| is the lattice constant, 
and taking into account that 

kri = 0, kr2 = — , krg = — , kr4 = 2kb, 



krg = — , kre = y , krr = 2kb, krg = ka = 3kb, (B2) 

ifo in k space becomes of the form 

^o = EEC'k,.^k4., (B3) 

o- k 

where Cl. ^ represents the eight component row vector {c\ k 0-5 4 k o-' •••; 4 k o-)' ^k.o- represents 
the eight component column vector obtained as the conjugate transpose of C*ko-' ^^^ ^k is 
a 8 X 8 Hermitian matrix whose transpose M^ is given by 
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FIG. 11: Exemplification of the bare band structure at ^2/^1 = 1.2,t3/ti = 1.3, t4/ti = 1.05,t']^/ti = 
0.1,^2/^1 = 0.2, eo/ti = — 0.9,ei/ti = 1.15, and ti > 0. The energy values are in ti units. Since 
E{k = ka) is even in k, only the k £ [0, it] is shown. The folding at k = it is a consequence of the 
structure of the cell. 
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FIG. 12: The bare band structure in the presence of the condition (jAllj) . One has ts/ti = 
1.3,t4/ti = l.05,t[/ti = 0.1,t'2/ti = 0.2,ei/ti = 1.15, h > 0, and t2,eo are given by (fATTI) . The 
energy values are in ti units. Since E{k = ka) is even in k, only the k G [0, ir] is shown. The 
folding at A; = TT is a consequence of the structure of the cell. 
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^^ikb/2 2t[ COS f tae" 



■ikb 



ei 



t^e 



Ikb 



ak 



-ikb 



2t'^cos^ tie'--" ei 



ofc 
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^ ikb/2 
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J. „ikb 
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tie 



-ikb/2 
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2t;cosf tge- 

t'l 
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^^^^kb/2 
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2t'2 cos f tie^^''/2 

t'l 

t; 

f^^-ikb/2 2t'2COsf 



eo 
toe 



toe' 



-jfc6 



eo 



t; 





t[ 

^ggifc6/2 2t; COS 



2t;cosf tge 



\ 



afc 



-ikb/2 



l^(,-ikb/2 2t;cos 



afc 



afc 



tse*^''/^ ^^g- 



ifcb 



t4e^'=^ 
ei 



;b4) 



2t\ cos 2 

The bare band structure is obtained from the secular equation of M^ as an eight order 
algebraic equation in A = En{k), n = 1,2, ..., 8 of the form 
1 



t2)2 



{{E' - Z' 



W 



|2\2 



where one has 

\F 
E 



ak 



Y' + 4Z[EY -Z{\W\' + \F\')]} = 0, (B5) 



ak 



^ = \V\^ + {\W\^ + 2X) cos' y , Y = 2{\W\^ + X) cos y , 

,2 .2\ „ /-/2 , ,2\ AJ-' U' J- J- \ „^„2 ^k 



eo(ei - t4) - ei(ti + tg) - 4ti(tiei - t3t4) cos — -, 
Z = 2[t'^{el - tl) - 2tf ei + t3t4t'i] cos — , X = d{a -b)+ AUt['^{a - 6 cos' — ), 

|Ty|2 = (|2 + 8t4t;'cCOs'y, 



IV'I 



a' 



h{h — 2a) cos' — , 



a 



t2{el - tl) + t4(tf - tl), b = ti(e? - tl) + 2t3(t;ei - t3t4), 



h {el - t^) - 2t; (t3ei - t;t4) , d = ti{ 



tiiei 



2t3t[ei, 



eo = eo - A, 



ei 



ei-A. 



(B6) 



An exemplification of the bare band structure is presented for the general case in FigJTTl 
while the bare band structure in the presence of the condition flAllI) is exemplified in FigIT2l 
As Figs.( fTT|T2l) show, fiat bands in the bare band structure are not present for non-zero values 
of all Ho parameters. This fact can be directly and analytically verified in the following way. 
The relation (IBSp can be written as 



ak 



ak 



As cos"" -^ + Aq cos" -^ + A4 cos^ 



ak 



A2 cos 



ak 



An = 0, 



(B7) 



2 " 2 2 2 

where the numerical prefactors A2m, f^t = 0, 1, ...,4 are dependent only on Hq parameters 
and the eigenvalue A. Since now in (IB7p the entire k dependence is concentrated in the 
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cos^™(A;a/2) terms, it results that flat band solutions (i.e. /c-independent expressions for A) 
can be obtained only if for a given solution, simultaneously for all A2m, ^2m = holds. 
But, since one has A^ = (16t| )^, this condition cannot be satisfied for non-zero value of the 
hopping matrix elements. Consequently, the band structure created by Hq does not contain 
fiat bands at all. One further notes that the A2m<8 prefactors can be as well easily calculated 
based on ( 1B5|B6p . and one finds that the denominator in f IBSP simplifies in all A2m- One 
obtains for example Aq = {{ej - t^)(eg -tj- tl) + (tf - tlf + 2[2tit-it\ei - eoei(tf + tl) - 
t2t4(tf - tl)]Y, or As = -256tf {t;^t'/e2/2 - [^'2^4(2^2^3 - t'^h) + t[^{2tit'^ + t2t3)]eit; - 
t'^%UeQ + t'l'eoei + t;^(tf + t^) + t'^htl^t'^h - t[h) + {t[\/A)[U{tl + 21'^^) + 2t3(4t;t'2 + 
2M^ + t2t3)-4tf(ti-t2)]}, etc. 

Appendix C: Example for W operators 

In this section one exemplifies solutions of ( TT7|) in the m = Q case (i.e. the number of 
cells in the test chain is six) inside the parameter space region specified by 

t[ > 0, t3 = 3t;, 6tf = U{ti - AQ. (CI) 

One further introduces the parameter A = {ti — 4^2) /(ti + 2t2)- 

In this region one has a B'^ solution holding the same coefficients in each cell, the solution 
being present at arbitrary Nc = m and arbitrary A. This operator will be denoted by the 
i = 1 index, and has the expression 

Afc-l 

-°l,o- ^ 2-^ ('^j+na+r5,o- ~ C^+na+r4,o- + '^j+(n+l)a+ri,cr ~ C^+(n+l)a+rfi,a>- l^^) 

n=0 

Given by its cell homogeneous nature, the -B^ from (IC2p cannot be translated or rotated in 
providing new, linearly independent 5^ operators. 

The following linearly independent solutions emerge for different A values from which 



32 



one presents below the A = \/3 case. For this parameter region, BI ^ has the coefficients 
_1 + V3 _ 1 + 73 _ 1 + ^3 _^ _ ^ 

2^2,1,1 — ^ ) 2:2,1,2 — ^ ) 3^2,1,3 — ? 2:2,1,4 — -L, 2:2,1,5 — "J-, 

2 2^^ 

1 + ^3 3 + VS 3 + ^3 



2^2,1,6 — 7^ — ) 2:2,1,7 — — —_ — ) 2:2,1,8 

2:2,2,1 = V 3, 0:2,2,2 = ) 2:2,2,3 = 7: ) 2:2,2,4 = — 7: — 5 2:2,2,5 



2 ' ^'^" 2g ' ^'^'^ 2g 

1 + VS _ 1 + Vs 3-Vs 3-Vs 

,2:2,2,3 — 7: — 

9 2^ 



/- 3 + V/3 

2^2,2,6 — ~V 0,0:2,2,7 — — ;:; ,2:2,2,8 — u, 

• -.-....^=-£-.~=^^'- -...=--.. 

-1 /^ -n - 3 + V3 

2^2,3,6 — J- — V 0,2:2 3 7 — U, 0:2,3,8 — ;:; , 

2^ 

3-V3 _l + v^ _l + y3 _ /- _ /- 

2^2,4,1 — ;:; ,2:2,4,2 — — ;:; ,2:2,4,3 — ,2:2,4,4 — —V 0,2:2,4,5 — Vo, 

2 2^^ 

_ 3-^3 _ 3 + ^3 _ 3 + ^3 

2^2,4,6 — 7: ,2:2,4,7 — 7: ,2:2,4,8 — ;:; , 

2 2g 2g 

_i + Vs _i + Vs _ 1 + Vs _i + Vs 

2^2,5,1 — "1,2:2,5,2 — ,2:2,5,3 — — ;:; ,2:2,5,4 — ;:; ,2:2,5,5 — 7: , 

g 2g 2 2 

-1 - 3 + V3 _ 

2^2,5,6 — -•-,2:2,5,7 — ;:; ,2:2,5,8 — u, 

2^ 

-n _1+V3 _ 1 + V3 _ _ 

2^2,6,1 — U, 2:2,6,2 — — ;:; ,2:2,6,3 — ;:; ,2:2,6,4 — u, 0:2,6,5 — u, 

2^ 2g 

2^2,6,6 = 0,X2,6,7 = 0,X2,6,8 = 7: , (C3) 

where g = 6t[/{ti + 2t2) is arbitrary. This solution can be translated five times by a ob- 
taining Bl ^, BI ^, ..., B^ ^, and the six operators BI ^, ..., B^ ^ can be rotated by vr obtaining 
-Bgg., ..., -B|3^. Consequently, at this moment, the number of linearly independent fi''' oper- 
ators with fixed spin at our disposal is 13. 
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The following linearly independent solution holding the index i = 14 has the prefactors 

2^14,1,1 = V3,a;i4i2 = ,3^14,1,3 = ,3^14,1,4 = 2 — V 3, X14 i 5 = —(2 — v3), 

9 9 

2^14,1,6 — ~V3,X1417 — — ,Xi4i8 — -, 

9 9 

/3 

3^14,2,1 = -2(1 - VS), a;i4,2,2 = , a;i4,2,3 = 0, Xi4,2,4 = 3 - 2a/3, Xu^2,5 = -(3 - 2^3), 

9 

1 

> 2^14,2,8 — 5 

9 9 



2^14,2,6 = 2(1 - \/3),a;i4,2,7 = 7, 



3^14,3,1 = -(3 -2\/3),Xu^3^2 = 0,Xi4,3,3 = — , a:i4,3,4 = 2(1 - y3),Xi4,3,5 = "2(1 - VS), 

9 

2 

.3:14,3,8 — — ) 
9 9 



2^14,3,6 — 3 — 2v3,Xi4 3 7 — , 

9 



2^14,4,1 — ~(2 — v3),Xi4 4 2 — — 53:1443 — — ,a;i4 4 4 — — v3,Xi4 4 5 — v3, 

9 9 

2^14,4,6 = 2 — v3,Xi4 4 7 = , 3:144 g = , 

9 9 

3^14,5,1 = 0,Xi4 5 2 = ,3:1453 = 0,a:i4 5 4 = —1, 3:145 5 = 1, 

9 

_ 1 _ 1 

3:14,5,6 — 0,a:i4 5 7 — — ,3:14 5 8 — — , 

9 9 

Vs 

3:14,6,1 = 1,3:14 6,2 = 0,Xi4 6,3 = ,3:14,6,4 = 0,Xi4 6,5 = 0, 

9 
1 2 

3:14,6,6 = "~1, 3:14,6,7 = -,3:14,6,8 = — • (C4) 

9 9 
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Similarly, another linearly independent BJ^ term £ = 15 is the following one 

_3-y3 _i-y3 _i _ rn r\ -or 

2^15.1,1 — 7: ) 3;i5 1 2 — ,3^15,1,3 — ; ^^is i 4 — —(2 — V 3), Xis 1 5 — 2 — V 3, 

2 2g g 

_ 3-V3 _ l-VS _l + V3 

2^15,1,6 — 7: )3;i5i7 — , 2;i5,l,8 — 7: ; 

2 2g 2g 

_„/^ _1 _ I + VS _5-3V3 _5-3V3 

3^15,2,1 — O — VO, Xi5^2,2 — 5 a;i5^2,3 — ^ ) 2:15^2,4 — ^ 5 3;i5_2,5 — ^ ) 

- rq /^^ - ^ + V^ _ 1 

2^15.2,6 — ~{'J ~ yo,Xi5 2,7 — 7: ; 2:15 2,8 — ", 

2t/ 9 

3^15,3,1 = 3 — V 3, Xi5^3^2 = 5 a;i5_3^3 = , a;i5^3^4 = —(2 — V3), Xi5^3^5 = 2 — V 3, 

9 ^9 

_3-v^ _ 1-V3 _ 1 _ ^ 

2^15,4,1 — 7: ,2:15,4,2 — — ;:; ,a;i5 4 3 — -,0:1544 — —i, 0:154 5 — i, 

2 2g g 

_3-V3 _1-V3 _1 + V3 

3^15,4,6 — 7: — ,2:1547 — — ,a:i5 4 8 — , 

2 2g 2g 

_1 _i + Vs _ I + VS _1 + V3 

3^15,5,1 — U,Xi5 5 2 — -,3:1553 — ,a;i5 5 4 — ,0:1555 — , 

g 2g 2 2 

-n - ^ + v^ _ 1 

2^15,5,6 — U,Xi5 5 7 — ,2:155.8 — , 

2^ 9 

-n _l + v^ _ 1-V3 _ _ 

2^15,6,1 — U,Xi5 6,2 — 7: ,a:i5 6,3 — ;:; ,2:15,6,4 — —-1-, 3^15,6,5 — -•-, 

2^ 2g 

..«^0,.,,„^-i,..,,^-i-^. (C5) 

9 ^9 

Both BI^^, and -8150- can be translated and rotated in providing new B^^ operators in the 
ground state wave vector, but from these not all are linearly independent. 

One further note that maintaining (JCip . solutions can be found also for other A values, 
as for example A = 0.75, 1, 2, 4, \/2, (v^ ± l)/2, etc. Furthermore, solutions are present even 
if (1C1|) does not hold, but in this case cell homogeneous 5"^ operators as presented in (1020 . 
are not present. 
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